Phase diagram of the t — t' — U chain at half filhng 
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We study the half filled Hubbard chain including next- nearest-neighbor hopping t' . The model has 
three phases: one insulating phase with dominant spin-density-wave correlations at large distances 
(SDWI), another phase with dominant spin-dimer correlations or dimerized insulator (DI), and a 
third one in which long distance correlations indicate singlet superconductivity (SS). The boundaries 
of the SDWI are accurately determined numerically through a crossing of excited energy levels 
equivalent to the jump in the spin Berry phase. The DI-SS boundary is studied using several 
indicators like correlation exponent Kp, Drude weight Dc, localization parameter zl and charge gap 
A,. 



INTRODUCTION 

Superconductivity in low dimensional systems has 
been a subject of interest in the last years. The physics of 
the negative-?/ Hubbard model, Q as well as that of the 
ladder Q suggest that the opening of a spin gap As 
is a key ingredient to superconductivity. In fact, a non 
zero As suggests the existence of bound pairs of electrons 
(like Cooper pairs), which if mobile can lead to a SS state. 
The indications of resonance valence bond superconduct- 
ing states in effective low energy models for 2D cuprates, 
like the generalized t — J 0,13 or the Hubbard model 
with correlated hopping, P are also consistent with this 
picture. Closely related with the latter system is the 
Hubbard model with an extra nearest-neighbor exchange 
J,|3] for which also strong signals of superconductivity 
have been found. 0, 

The half-filled Hubbard model with additional next- 
nearest-neighbor hopping t' in one dimension is appeal- 
ing because, for weak U , one expects the existence of 
a metallic phase when t' > 0.5i,[3)E| and, from pre- 
vious studies of the frustrated Heisenberg chain. [l^ 
one expects the opening of a spin gap for large U and 
t' /t > with Rc = 0.241 167. In addition, this 

model might be valid as an effective model for several one 
dimensional cuprates. ,14.. .15.] The corresponding Hamil- 
tonian is: 

■i(T i(T 

+U^n,^nii. (1) 



sketched an approximate phase diagra m, lid confirmed 
the SS character in the metallic phase, [l7| and deter- 
mined some points where appro ximately the charge tran- 
sition SS-DI takes place. [l7llll| 

In this work we present further evidence of the domi- 
nance of SS correlations in the metallic phase by calcu- 
lating the correlation exponent Kp. We also obtain more 
information on the metal-insulator transition using nu- 
merical calculations of the Drude weight Dc, localization 
indicator 2;^ [l9ll2oll2]| and charge gap Ac. Furthermore, 
we determine the opening of the spin g ap As u sing the 
crossing of singlet and triplet excitations, [lj,|22| which is 
equivalent to the jump of the spin Berry phase. ^23^ This 
method has been shown before to lead to very accurate 
results for this transition, even if the size L of the s yste ms 
considered is less than IS.dlHllllilllllilliili^ As 
an example, the resulting opening of the sp in g ap for the 
Hubbard model with correlated hopping |3, 123 coincides 
with field theory results in the weak coupling limit ■fH.lBOl 
and with exact results when the number of doubly occu- 
pied sites is conserved, js^ Moreover, this method allowed 
us to determine the complete boundary of the SDWI 
phase. 

In Section II we present the results for the spin transi- 
tion and briefly describe the method used. In Section HI, 
we discuss our results for the correlation exponent Kp. 
In Section IV we explain briefly the meaning of the local- 
ization indicator z^, present results for it, for the Drude 
weight Dc, and for the charge gap Ac, and use this in- 
formation to draw a tentative a phase diagram. Section 
V contains a summary and discussion. 



Using weak coupling renormalization group, Fabrizio 
obtained that a spin gap opens for t' > 0.5t, and that 
for As ^ SS correlations dominate in the metallic 
phase while spin dimer correlations are the largest at 
long distances in the insulating phase, [loj Later stud- 
ies using density matrix renormalization group (DMRG) 



THE SPIN TRANSITION 

This transition can be determined by a topological 
transition (a jump in the spin Berry phase), which is 
equivalent to a suitable crossing of excited energy lev- 
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FIG. 1: Topological sectors of the spin Berry phase in the f', U plane 
(t = 1). As = for 7s = 7r and As / for 7, = 0. 
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FIG. 2: Same as Fig. 1 in the 1/U, t' plane. The open symbol 
corresponds to the known result t' — \/Rc — 0.4911 for U — > cjo. 



els. The basic idea of this crossing is that, in confor- 
mal invariant systems (corresponding to fixed points in 
weak coupHng bosonization theory plus renormalization 
group), the lowest excitation energy determines the dom- 
inant correlations at large distances. [l^E3 | On the other 
hand, the topological transitions denote rearrangements 
of charge or spin which have a very transparent mean- 
ing in the strong coupling limit, [2J, |23, |29j and are ex- 
pected to be valid by continuity also for weak coupling. 
For SU(2) invariant systems both criteria coincide for the 
spin transition. 23] 

The spin Berry phase 7^ is the phase captured by the 
ground state using twisted boundary conditions -e'* for 
spin up and e~'* for spin down- in the cycle < $ < 
27r.|2^] In gauge invariant form, it can be numerically 
calculated as: 

7s = - lim {Im[ln(n;i-2^.g(a>„-ci>,)|g($,+i,-$,+i)) 
(g($^_i,-<i>jv-i)|g(27r)))]}, (2) 

where is the ground state of -ff and 

\g{2n)) = exp[i^^pX"^T " ^jl)M^M- The Hamil- 
tonian H differs from H in that the hopping term has 
the form -tY.^M+laCi.e'*"''' + H.c). 

The Berry phase 7s can take only two values: or tt 
(modulo 27r). The jump between these values is related to 
the change in the difference between electric polarization 
from spin up and down by A7S = (27r/e)A(P| — P[ V|2^ 
The position of the jump of 7s coincides with the singlet- 
triplet crossing of the energy levels: the crossing of a 
singlet even under inversion with an odd triplet, with 
K = 0. The boundary conditions are periodic (antiperi- 
odic) for system lengths L — 'in (4n -I- 2), with n inte- 
ger (they are opposite to the so called closed shell ones). 
The same crossing is the one which determines the open- 
ing of the spin gap (the Kosterlitz-Thouless transition 
in the spin sector) according to conformal field theory, 



after bosonizing the problem and using renormalization 
group. 22] This theory predicts a l/Li^ scaling of the tran- 
sition point. 

We have calculated the transition for L = 8, 12 and 
16. The restriction to L values multiple of four is due 
to the fact that other system sizes lead to frustration of 
the spin-density waves (or classical antiferromagnetic or- 
dering) expected for large U . In particular, for i = the 
even and odd sites are decoupled and, to avoid frustration 
of the AF classical order in each subsystem, L/2 should 
be even. We have fitted the resulting three points with 
a straight line a + hjl? in Xjl? . The transition point in 
the thermodynamic limit is determined by a and its error 
gives an indication of the fit quality. 

In Fig. ^ we show the results for the spin transition 
for [/ < 10 (we have chosen t = 1 as the unit of energy). 
For t' < 0.5 we obtain that the transition takes place at 
U — for all L. For positives values oi U, 7s = tt and 
As = 0, while for negative U, 7s = and As ^ 0, as 
in the ordinary Hubbard model. |23 At the point U — 0, 
t' — 0.5, the spin transition line deviates 90°, and for 
positive U its shape is like a nose, with a maximum value 
of t' = 0.564 ± 0.001 taking place near U = 4.5. The 
errors of the points shown in Fig. Qare less than 0.003. 

Results for larger values of U are shown in Fig. |21 For 
t, t' ^ U, to leading order in 1/U the model is equivalent 
to the Heisenberg chain including next nearest-neighbor 
exchange: 

if = ^(JS,.S,+i+ J'S,-S,+2) (3) 

i 

with J = At'^/U, J' = W^/U. Previous studies of this 
model using the method of crossings of excitation levels 
plus a detailed knowledge of the size dependence have 
established the critical value Rc = {J' / J)c — 0.241167 at 
the transition. Then, for U 00, As should open at 
i'c — VRc- Our results are consistent with this fact. A fit 
of our transition points at U = 15, 20, 30, and 50 with 
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the parabola t'^ — \fWc = a + hjU + cjU"^ gives c = 3.6 
and a = & = within our numerical error for t'^ (less than 
5 X lO^**). The fact that 6 = is also consistent with the 
strong coupling expansion: h is proportional to the first 
correction to Eq. which is of order t^t' /U"^ . However, 
the electron-hole transformation cj^. ^ Cia changes the 
sign of both hoppings but should leave invariant the spin 
gap at half filling. Then, all the corrections to H of order 
1/^/2"- with n integer should either vanish or be irrelevant 
for the spin gap since they change sign under electron- 
hole symmetry. 



CORRELATION EXPONENT IN THE METALLIC 
PHASE 

In the metallic phase, the analysis based on bosoniza- 
tion and renormalization group 10] shows that there is 
only one gapless charge model and the spin modes are 
gapped (A3 7^ 0). This analysis is valid at weak coupling 
and requires that t' > 0.5. [ll| In this case, there are four 
Fermi points at wave vectors ±fci?i and ±fci?2- The decay 
of the different correlation functions at large distances is 
determined by the correlation exponent Kp. Depending 
on the value of this exponent, the slowest asymptotic de- 
cay corresponds either to spin dimer-dimer correlation 
functions or to superconducting ones. At half filling the 
former have the following decay: 



XDw{.x) = 



where 



{ODw{x)ODWm 
C0S{'KX) 



cos[2(fcir2 — kpijx] 



(4) 



Odw{x) = S+{x)S-{x + a) + S+{x + a)S-{x). (5) 
For the dominant superconducting correlations one has 



Xscix) = (A(x)At(O)) 
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(6) 
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FIG. 3: Correlation exponent Kpi defined by Eqs. as a function 
of 1/L, for several values of U and t — 2.5. Triangles Kpi, circles Kp2, 
squares iCps. 



the charge velocity, k is the compressibility and Dc is the 
Drude weight. In finite systems these quantities can be 
calculated as: 



E{K = 2tt/L) - E{K = 0) = 



27r 



(9) 



where 



^ E{N + 2)+EiN-2)-2E{N) ^ ^^0^ 



Mx) = ^Pknfix) ■■ i^-pkfUi^) 

-i'pkF2T{x) ■■ 'ip-pkF2lix)- 



(7) 



An analysis similar to that of Schulzj31| extended to 
this case with 4 Fermi points leads to the following rela- 
tions involving the correlation exponent Kp: 



1 



4:K. 



'2vpKp2 



pi 



(8) 



Here, the subscript 1 or 2 in Kp is due to the fact that, 
in finite systems, when Kp is calculated using these ex- 
pressions the results differ. In the above equation Vp is 



Ld^Eo{L,'^>) 
2 



a$2 



(11) 



where E{K) is the energy of the lowest singlet with total 
wave vector K, E{L, N) is the ground state energy for 
length L and N particles, and Eo{L,^) is E{L,L) in 
presence of a flux $/i/(27re) threading the ring. Unless 
otherwise stated, the boundary conditions correspond to 
the "closed shell" ones (those which lead to the minimum 
energy): <& = (tt) for N = 2n with N odd (even). 

The above equations lead to two independent ways of 
calculating Kp. A third possibility uses the geometric 
mean Kp^ = ^Kp\Kpi. In Kp-^ the charge velocity does 
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FIG. 4: Localization indicator zi^ and Drude weight as a function 
of U for several system sizes h and t' — 2.5. 



not enter, and usually it leads to a more accurate result. 
The consistency of the three results in the thermody- 
namic limit is a test on the validity of the Luttinger liq- 
uid description, with gapless charge excitations. In fact, 
when the Kpi extrapolated to the thermodynamic limit 
start to deviate from each other one has a strong indica- 
tion of a metal- insulator transition. 32] For example, in 
the model with correlated hopping, the point where the 
Kpi start to deviate agrees with the jump in the charge 
Berry phase which corresponds to this transition. 27]. i32ij 

In Fig. 13 we show the evolution of Kpi as \J increases, 
for t' — 2.5. Clearly, for C/ = 10 the extrapolated Kp are 
inconsistent and the system is in the insulating phase. 
The metal- insulator transition seems to be near \J jt = 6, 
where Kp ^ 1/2. For lower values of [/, Kp increases, 
which indicates that superconducting correlations are the 
dominant ones inside the metallic phase. This is in agree- 
ment with previous results, lid, llj 



THE CHARGE TRANSITION 

In absence of disorder, the Drude weight Z?c defined in 
the previous section is an indicator of the metal-insulator 
transition. In particular, in the thermodynamic limit it 
should go to zero in the insulating phase, [l^l In recent 
years, Resta proposed to use as an alternative a localiza- 
tion indicator zr,.|,19.J However, a simple symmetry argu- 
ment shows that this expression gives an incorrect result 
for periodic systems with non integer number of parti- 
cles per unit cell One of us and Ortiz have gener- 
alized this concept for correlated systems with rational 
n,|23,0| following a corresponding generalization for the 
charge Berry phase. In the present case, with n = 1 
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FIG. 5: Same as Fig. 4 for t' = 5. 



one has 



ZL = 



(12) 



and zl (1) for a metal (insulator) in the thermody- 
namic limit. Note that the operator entering Eq. H12|l is 
a displacement operator in momentum space. Then, for 
an extended system with a well defined Fermi surface this 
operator shifts the Fermi surface and the scalar product 
vanishes (z^ = 0). Instead, for an uncorrelated band 
insulator the operator reproduces the same state and 
zl — I- Other cases with either disorder or correlations 
were also discussed before, [iol El| One advantage of zj, 
over Dc is that it can be calculated with DMRG.|2ll.E7l| 
although we are not considering this here. 

In Fig. 01 we show the results for Dc and as a func- 
tion of U, for t' = 2.5 and L — 8,12. Both quantities 
vary smoothly with U. While the results point towards 
an insulating behavior for U > At' = 10 (in particular 
zl > 0.5 and increases with system size) and a metal- 
lic behavior for small U (z^ and decreases with 
L for U < 1.75t'), it is not possible to determine ac- 
curately the critical value of U for which the transition 
takes place. For large system sizes there is a level cross- 
ing in the ground state, and while for one of the states 
zl increases with L, for the other one zl decreases with 
size. This suggests to associate the crossing point with 
the metal-insulator transition, which gives Uc — l.75t'. 

As showed in Fig. |5l when t' is increased to i' = 5 the 
values of Dc increase and those of zl decrease, suggesting 
an enlargement of the metallic region. However, from the 
dependence on L of the corresponding quantities no safe 
conclusions can be drawn about the critical value of U 
at which the transition takes place. For the same ratio 
t' = 5 we have also calculated numerically the charge 
gap, including DMRG results up to 20 sites, using the 
expression: 

Ac = E{L, N + l) + E{L, N -1)- 2E{L, N). (13) 
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FIG. 6; Charge gap as a function of 1/L for i' — 5 and different values 
of U. From bottom to top: 2, 6, 8, 8.5, 9, 10, and 20. FIG. 7: Phase diagram of the t — t' - U model. The star and triangle 

were taken from Refs. Ilgll and 1171 respectively. 



The results are shown in Fig. ?? together with a hnear 
extrapolation to the thermodynamic limit. Although a 
straight line in 1/L fits well the results for all values 
of U, the extrapolated Ac becomes negative for U < 
6, indicating again that finite size effects are important 
(since Ac should be either zero or positive). Nevertheless, 
for [/ > 10 it seems clear that the gap is open. 

Unfortunately, the charge Berry phase does not jump 
at the metal-insulator transition, in contrast to the case 
of the Hubbard model with correlated hopping. 0, E^l 
The study of this transition using Kpi, Dc, and Ac as 
explained above is not enough to lead to a well defined 
phase boundary. With this information we can draw a 
tentative phase diagram, presented in Fig. [7| with large 
error bars for the metal-insulator DI-SS boundary. We 
also include some points obtained by other authors us- 
ing larger systems. 0, 0| The boundary of the SDWI 
determined in section II is also plotted. 



SUMMARY AND DISCUSSION 

We have studied the phase diagram of the t — t' — U 
model at half filling, using different quantities calculated 
mainly by numerical diagonalization of finite rings. Un- 
fortunately, this model has large finite size effects in the 
regime in which there are four Fermi points in the non 
interacting case {t' > 0.5t). This can be checked by calcu- 
lating correlation functions or other quantities as a func- 
tion of size L in the non interacting limit, either for a ring 
or with open boundary conditions. Kuroki et al. argue 
that one should choose system sizes for which the level 
structure of the different bands in the non interacting 
case is similar. [T^ 

The above mentioned difficulties prevented us to ob- 
tain quantitatively valid results for the metal-insulator 
transition. Nevertheless, our calculations of Kp con- 
firm the dominance of superconducting correlations in 



the metallic phase found earlier by analytical [lOj and 
numericalflTj methods. Our tentative phase diagram also 
agrees with previous studies. (l^ilTnl^] For the spin tran- 
sition, on the contrary, our results show the 1/L^ scaling 
expected from conformal field theory^, |23| and, there- 
fore, are reliable. Moreover, they converge with high pre- 
cision to the known result in the strong coupling limit 
up to terms of order l/U"^. The jump of the spin Berry 
phase and the equivalent level crossing method have been 
successful in determining the openin g o f the spin gap in 
different SU(2) invariant models, 0, 111 [H IH El El 
and its validity has been confirmed by other analytical 
and numerical methods. 0, 113 In the present case, we 
became aware of recent DMRG results for this transition 
based on extrapolations of the dimer order parameter. |35j 
These results are very similar to ours except for a shift to 
larger t' by about ^ 0.7i for U ~ bt. We believe that this 
might be due to difficulties in detecting an exponentially 
small dimerization order parameter or gap using DMRG. 

An interesting question is what is the origin of the 
superconducting correlations in this system, where the 
only interaction is repulsive. An intuitive picture can 
be built starting from the strong-coupling limit C/ — > cx), 
where the model is described by the spin Hamiltonian Eq. 
(PI, whose properties are well known. 12] In the dimer- 
ized phase, at the Majumdar-Ghosh point J' = 0.5 J the 
ground state consists of nearest-neighbor spin singlets. 
These electron pairs are static and the system is insulat- 
ing. However, when U is reduced these pairs can overlap 
and become mobile, behaving similarly to bosons with 
short range repulsion, for which dominance of singlet su- 
perconducting correlations at large distance can be ex- 
pected. 
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